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Abstract 

We discuss in detail the calculation of the leading order QCD 
corrections to the Effective Hamiltonian which governs 6^57 and 
b ^ s g transitions in three different regularization schemes (HV, NDR 
and DRED). We show that intermediate stages of the calculation do 
depend on the regularization, but the same scheme independent co- 
efficients can be obtained in all the considered schemes. A detailed 
discussion of the results already present in the literature is also given. 
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1 Introduction 



Leading order (LO) QCD corrections to the Effective Hamiltonian governing 
the radiative decays of the B meson have been calculated by many authors in 
the last five years Ii^ fact they turn out to be large and important for 

the B phenomenology. In spite of this effort, different results are still present 
in the recent literature on this subject, even if the residual differences are 
numerically small and have essentially no phenomenological relevance. The 
aim of this paper is to help clarifying the origin of these differences by giving 
full details on our calculations in three different regularization schemes. 

Let us briefly recall how the calculation of LO corrections to the Effective 
Hamiltonian for radiative B decays has been developed in the recent years. 
The original calculations, refs. ^ and p|, were performed in the Naive Di- 
mensional Regularization scheme (NDR) and in the Dimensional REDuction 
scheme (DRED [H) respectively and used a reduced set of operators for the 
Effective Hamiltonian. They disagree on the results for the anomalous di- 
mension matrix, that was believed to be regularization scheme independent 
at the leading order. 

Later NDR result was confirmed [Q] and this led the authors of ref. 
to cast doubts on the reliability of the DRED scheme. Two years ago a first 
paper appeared where the complete LO correction was calculated in the 
NDR scheme using the full set of operators, then other similar calculations 
followed [3-0. Refs. and confirm the result of ref . |^ for the anoma- 
lous dimension matrix in the "reduced" basis. However they disagree on 
some new matrix elements in the "full" basis. 

Recently the question of the scheme independence of the Effective Hamil- 
tonian for B meson radiative decays has been clarified [0. The calculation 
of the complete LO corrections has been performed in NDR and in the 't 
Hooft-Veltman scheme (HV [T^) and it has been shown that the final results 
for the Wilson coefficients are regularization scheme independent provided 
one takes properly into account the scheme dependence of the one and two 
loop Feynman diagrams. Calculations which use the reduced set of operators 
such as those of refs. [0, |], Q have been demonstrated to give results which 
depend on the regularization scheme. Unfortunately the new NDR anoma- 
lous dimension matrix obtained in ref. differs from both the results of 
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refs. 

In this paper we provide full details on the calculation presented in ref. 
[0], including tables with the pole coefficients for all the diagrams. We also 
present the calculation in the DRED scheme. We show how the same reg- 
ularization independent results can be obtained in DRED as well as in the 



other two schemes, thus extending the result of ref. [jTT|. We also discuss 
why the conclusions of ref. on the reliability of the DRED scheme are 
wrong. Finally we compare our results with the most recent ones ^ and 
comment on the differences among the results, reported in refs. [pT| , |12|. 

The paper is organized as follows. In Sec. ^ we recall the Effective 
Hamiltonian for 6 — s 7 decays and the renormalization group equations 
(RGE) governing the evolution of the Wilson coefficients. The full operator 
basis is given and compared with the reduced one. The initial conditions for 
the coefficients are also given. In Sec. | the explicit solution of the RGE, 
showing the scheme independence of the coefficients, is given. In Sec. ^ we 
summarize the definitions of NDR, HV and DRED regularization schemes, 
discussing how they are implemented in our calculations. In Sec. ^ the 
contributions of one loop diagrams to anomalous dimension, operator matrix 
elements and counter-terms are considered and diagram by diagram results 
are given. The same is done for the contributions of the two loop diagrams 
to the anomalous dimension. We discuss the relations among one and two 
loop diagrams induced by the scheme independence of the final result and 
briefly report other checks done on our calculation. In Sec. |^ the anomalous 
dimension matrices in the different schemes are given. Finally in Sec. we 
critically compare our results with the other ones present in the literature. 

2 Effective Hamiltonian: Initial Conditions 
and Evolution for the Coefficients 

The effective Hamiltonian for 6—^37 [b —>■ s g) transitions is given by 

Heff = -Vt,V*^j2Q,{f,)Q{f,) ~ g^(/i)C(/i) (1) 
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where Vij are the elements of the CKMpTSI, quark mixing matrix. We use 
the following operator basis Q 

Ql = iSaCi3)(V-A)iCi3ba)(V-A) 

Q2 = iSaCa)(V-A)iCi3bf3)(^V-A) 

Q3,5 = {Saba){V-A) ^ {Wp)(VtA) 

q=u,d,s,--- 

Q4,6 = {Sabi3){V-A) i.Wa)(VTA) 



=u,d,! 



— -^^'^bSa(^'(v+A)^tf3bl3Gf,y. (2) 

Here {V =F ^) indicate the chiral structure, a and j3 are colour indices, m;, 
is the b quark mass, = — | is the electric charge of the down- type quarks 
and g (e) is the strong (electro-magnetic) coupling. The colour matrices 
normalization is Tr{t^t^) = 5^^ /2. 

The choice of the operator basis deserves some comments. In fact the 
basis in eq. (^) is obtained by using the equations of motion of the external 
fields. While this procedure has been criticized in the past, a recent paper 

shows that it can be safely used. This basis is often called the "complete" 
basis, to be compared with the "reduced" one of refs. [||, ||, Using the 
"reduced" basis, one neglects the contributions of the operators O3, . . . , Og. 
Also retaining the penguin operator| Op = S'^^D^t'^G^J) as in ref. 0, one 
actually does not fully consider the contribution of O3, . . . , Og. Unfortu- 
nately just the insertion of O5, Oq in the one loop diagrams gives a scheme 
dependent contribution, which does not vanish in the NDR scheme. Hence 
NDR calculations using the reduced basis give a regularization dependent 
result, see ref. [0 and Sees. ^, ^, below. 

The coefficients C(/i) of eq. (|l]) obey the renormalization group equations 

d _2:!M)(7(t,«^(t)) = o, (3) 




where t = ln(M^//i^) and as = g"^ /^t^- The factor of 2 in eq. normalizes 



the anomalous dimension matrix as in ref. ||T5[. 7 includes the contribution 



^In the "complete" basis this operator can be ehminated via the equations of motion. 
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due to the renormalization of m?,, the gluon field and the strong couphng 
constant see e.g. refs. 0, |18| (see also Sec. |^). 

The initial conditions for the coefficients can be easily found by matching 
the effective theory with the "full" theory at the scale Mw- They are given 



by [O 



Ci{Mw) = 
C2{Mw) = 1 
C3{Mw),...,C,{Mw) = 



, 3x^-2x^^ 8x^ + 5x'^ -7x 



where x = ml/M^r. 



2(1 4(1 

x^ — — 2x 
2{l-xY 4(1 -x) 



C^{Mw) = —lnx + —^^ (4) 



Regular izat ion Scheme Independence of the 
Effective Hamiltonian 



In this section the explicit solution of eq. (|^) for the coefficients C(/i) is given 
and the regularization scheme independence of the Effective Hamiltonian is 
discussed, following ref. 0. The solution in eq. (^, as well as the relation 
among matrices in different schemes in eq. (0), rely on the peculiar structure 
of the anomalous dimension matrix, which is 
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as 
Air 



7r P7 
0" 777 







(5) 



V U 787 788 / 



This "almost" triangular form is obtained because the magnetic operators Qj, 
Qs do not mix with the 4-fermion operators Qi, . . . ,Qq. Then it is convenient 
to introduce the reduced 6x6 matrix 7^ which mixes Qi, . . . ,Qq among 
themselves and two 6-component column vectors = (717, 727, . . . , 757) and 
Ps = (7185 728, • • • 5 768), which account for the two-loop mixing of the 4- 
fermion operators with the magnetic ones. 
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A very peculiar characteristic of this calculation is that at the LO of 
the perturbative expansion both one-loop and two-loop diagrams contribute. 
The mixing of 4-fermion operators among themselves (7^) as well as the 
mixing of magnetic- type operators among themselves (777, 737, 788) are one- 
loop effects (at LO), while the mixing between 4-fermion and magnetic-type 
operators {P^, Ps) is a two- loop effect (always at LO). 

In terms of these quantities the RGE are given by 

2l,'^^CM = ^(/37.(7,(/x)+777C7(/x)+787^^8(/^ 

2/.^^C8(/x) = ^ {Ps ■ Cr(ii) + IssCs(l^)) , (6) 

where Cr{ij) = (Ci(//), . . . , C6(/i)), cts = cts{ij) and fji^d/dfi^ = fx^d/dfj? + 
P{as)d / das- Diagonalizing the submatrix which mixes the magnetic opera- 
tors, one obtains 

2/.^^a(/.) = ^^^-CM 

^ ^788'y8(/^), (7) 



where 



= C7(//)+a7-aM+ CM 

111 - 788 

vM = C8(/x) + a8-a(/x) (8) 



with 



«7 = (777! - 7r) 



Pi + ^^/38 

777 - 788 



0L% 



(7881 -7r) Vs- (9) 



5 



The rhs of eqs. (|^) involves only one loop quantities, so that the solu- 
tions Cr{[j), vji^fi) and vs{fi) are independent of the regularization scheme. 
However the expression of both f7(/i) and Vsi^), eq. (||), contains both 
scheme-dependent and -independent quantities. In particular C.r(/i) is the 
LO solution of the RGE for the 4-fermion operators, which is known to be 
regularization scheme independent. Indeed, Cr{n) is known up to next-to- 
leading order terms |]T3|, , which are not to be considered in our LO calcu- 
lation. On the contrary jSj and /^g, hence aj, cJg and C-r^fi), Csi/J^), come from 
two loop diagrams (see Sec. ^ below) and they do depend on regularization 
scheme, even if they are LO results. Precisely this point, that was missed in 
refs. is responsible for the difference between previous results in NDR 

and DRED schemes, as stressed in ref. 0. 

Moreover, already noted in ref. the operators Qs, have non zero 
matrix elements between the b and 57 {s g) states, through the penguin 
diagrams in fig. with massive b loop propagators. These one loop matrix 
elements and the tree level matrix elements of the magnetic operators are of 
the same order in as, so that the matrix elements of the Effective Hamiltonian 
can be written as 

<s-f\Heff\b> = Crijj,) < s-f\Q7{iJ,)\b > +C5{iJ,) < s-f\Q5{^i)\b > 
+C6{n) < sj\Qe{fi)\b > 
= Crifi) < s-^\Q7{fi)\b > 



<sg\Heff\b> = Csifi) < sg\Qs{n)\b > +C5{fi) < sg\Q5ifj,)\b > 

= Csifi)<sg\Qsifi)\b>, (10) 

where the coefficients C-j^jj), CsifJ^) are defined as 

Crifi) = CM + Zr ■ Crifi) 

Csif^) = C8{ii) + Zs-Cr{ix). (11) 



Eqs. (p!0| ) and (|TT]) can be seen as a finite renormalization of the operators 
such that the matrix elements of and Qq vanish. The vectors and Zg 
can be considered as the effect of a mixing of order a° among Qs, Qe and the 
magnetic operators. They are calculated from the finite part of the penguin 
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diagrams in fig. (j^), thus they depend on the regularization scheme. We find 
that they vanish in HV and DRED, while in NDR they are given by 



2^NDR 



(0,0,0,0,2,2A^) 
(0,0,0,0,2,0). 



(12) 



We now show how the scheme independence of the Effective Hamiltonian 
is recovered. In terms of the operators renormahzed at the scale /x, the 
Effective Hamiltonian is given by 

Q^/i) ■ a(/i) 

-1 / - 787 



H. 



eff 



+ 



787 



(7, - 7771) P7 + 



777 - 788 



a(/i) 



^^8 



(/i) + (7r - 7881) ^8 ■ CrifJ') 



788 - 777 
+ I^Vsifl)+ (7r- 788i) Ps + Zs 



787 



[ 788 - 777 

+ [vsifJ^) + UJs ■ Cr{f^)}Qs{fJ^), 



(13) 



where 



LJj = - 7771) + 

^^8 = {%■ - 7881) (38 + Zs. 



787 



777 - 788 



-Ps + Zj 



787 



777 - 788 



(14) 



Thus the Effective Hamiltonian is scheme independent provided that 007 and 
Us do not change with the scheme. We will prove in the following that this 
is indeed the case. 

Now we discuss the change in the anomalous dimension matrix induced 
by a change of the regularization scheme. The anomalous dimension matrix 
in a given scheme "a" is defined as 



r = 2(z'^(/i))~V2^z'^(/x), 



(15) 
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where Z°- is the matrix of the renormahzation constants which gives the 
renormahzed operators in terms of the bare ones 



Q(/i) = (Z'^) ^Qb. (16) 

Let us consider different renormahzation schemes. The operators renor- 
mahzed through the usual MS subtraction change from one scheme to an- 
other. In order to have the same renormahzed operators in aU the schemes, 
one can define them using the MS procedure in a given scheme, then adopts 
a suitable non-minimal subtraction in the other ones. So the MS renormal- 
ization constants in two different schemes "a" and "6" can be related through 
the equation 

= Z^r, (17) 

where the matrix f accounts for the change of the subtraction procedure in 
the scheme "6" necessary to define operators renormahzed as in the scheme 
"a" . In our case the matrix f is expressed in terms of the vectors Z introduced 
in eq. ([TT|), as discussed in the following. 

The relation between the anomalous dimension matrices in the schemes 
"a" and "5" is easily obtained from eqs. ([T5|)-(p^ 



This relation was already found in ref. |2y] and applied to the next-to-leading 



calculation of the AS* = 1 Effective Hamiltonian. In that case f deviated from 
1 by terms of order a^. Here we are not interested in such terms, but now, as 
a consequence of eqs. ([T0| ) -([T2| ) , f differs from the identity already at order 
However, due to the peculiar form of the anomalous dimension matrix 
discussed above, the matrix f is simply 

(19) 

where 16,2 are 6x6 and 2x2 identity matrices and AZ is a 6 x 2 matrix. 
Imposing the condition that the renormahzed operators coincide in the two 
schemes, AZ is defined as the difference of the vectors Z of eq. (0) calculated 
in the two different regularization schemes, i.e. AZ = {Zj — Zj, Z^ — Zg). 
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At the order we are interested in, eq. (|T^) gives 



A7 = 7" 



7 



[AZ, 7''] - AZ7''AZ. 



(20) 



Contrary to the aheady mentioned AS" = 1 case, eq. ( PD| ) does not define a 
scheme independent combination of the one and two loop matrices, namely 
it can not be written just in terms of quantities which are either calculated as 
differences between the two schemes or scheme independent by themselves. 
However a scheme independent combination actually exists due to the struc- 
ture of the matrices. In fact, applying eq. (^) to anomalous dimension 
matrices of the form of eq. @ , we find 



(A/3r), 
(A/38), 



A7,7 

A7j-8 



(777! - 7r) AZ7 + 787AZg 

(7881 - 7r) AZ8 



(21) 



where j = 1, . . . , 6. This clearly implies that the combinations 



/?7 - (777I - It) Zj - 7872'8 
/38 - (7881 - 7r) ^8 



(22) 



are scheme independent. Using eqs. (|2T|) , one can easily check that the com- 
binations tJy, tJs in eq. give Aci;7^8 = 0, i.e. they are scheme independent 
too. In turn this implies that, as expected, the Effective Hamiltonian, eq. 
(|13D, is independent of the regularization scheme. 

Eqs. ( pTD can be used to relate classes of diagrams calculated in different 
schemes. Thus it can be a useful check of the calculation, see Sec. |^ below. 



4 Dimensional Regularization Schemes 

In this section we briefly recall the deflnitions of the regularization procedures 
we have used for our calculations. Other nice discussions can be found in 
refs. pOl ^ and references therein, while for more details the reader should 
refer to the literature on the speciflc subject. 
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The regularization procedures we have used, NDR, HV and DRED, are 



all based on the dimensional regularization of the loop integrals [|T0], which 
consists in performing the integration over the loop momenta in D dimen- 
sions, thus turning the divergences into regularized 1/(4 — D) terms. The ex- 
tension of the Dirac algebra to non-integer dimensions presents no difficulty 
except for the properties of traces involving 75 and in general expressions 
containing the completely antisymmetric tensor e^^p^, which does not have 
any meaningful extension in D dimensions. Actually the three schemes differ 
from each other in the way they treat 75. Another relevant point is that the 
relation 

7^,7u7p = l^idup - luQpp + lpg^,u - ^1''lb€-^,vpa (23) 

which, in 4 dimensions, projects the product of three 7 matrices on the basis 
does not hold in D dimensions. This implies that more complicated tensor 
structures, which eventually vanish in 4 dimensions, appear in the regularized 
theory and must be considered in the renormalization procedure. 

Let us now separately discuss how the three schemes are implemented. 



4.1 Naive Dimensional Regularization 

In NDR all the Lorentz indices appearing in the regularized theory are as- 
sumed to be in D dimensions. The Dirac algebra is identical to the 4- 
dimensional one, including the properties of 75, once the 4-dimensional metric 
tensor is replaced by the D- dimensional one. The definitions we are inter- 
ested in are simply 

{7^,75} = 0. (24) 

Obviously, since the completely antisymmetric tensor is not defined in D- 
dimensions, traces involving odd number of 75 are ill-defined too. However, 
to our knowledge, once the problem is properly fixed, there is no calculation 
which shows a failure of this scheme. Our way to fix the 75 problem is the 
one implemented by Schoonschip |^2|, which defines traces of odd number of 



75 in NDR in the same way they are defined in the HV scheme, see below. 
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4.2 't Hooft-Veltman Regularization 



In HV [|T0|, Lorentz indices in D, 4 and {D — 4) dimensions are introduced, 
together with the corresponding metric tensors g^u, g^i, and g^^. All the 7 
matrices are taken in D dimension, then indices are split in 4 and {D — 4) 
components, according to the rules 

-g^.ug"' = 0. (25) 

These rules define also the extended Dirac algebra, once one notes that 7 
matrices in 4 (7'') and {D — 4) (7'') dimensions can be written in terms 
of D dimensional matrices as g^'^'ju and g'^^'~fu respectively and the usual 
commutation relations among the D-dimensional Dirac matrices in terms of 
the D-dimensional metric tensor gf^^ are assumed. 

7 matrices in D dimensions do not have definite commutation relation 
with 75. In fact the following relations hold in HV 

{%,75} = , [%,75]=0. (26) 

This is equivalent to define 75 as the product z7o7i7273- This way of treating 
75 in D dimensions is the only known one which does not give rise to algebraic 
inconsistencies or introduce ill-defined quantities. On the other hand, due to 
the splitting of indices, the HV scheme is the most difficult to handle among 
the three considered here as far as the algebraic manipulation problems are 
concerned. 

Finally we mention that the chiral vertices in D dimensions can be defined 
in different ways, all having the same limit when D tends to 4 dimensions. 
We use the symmetrized form 

^iV±Ah,{V±A)=%{V±A). (27) 

In this way the bare vertices preserve the chirality of the external fields also 
in D dimensions. 
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4.3 Dimensional Reduction 



In DRED PI the Dirac matrices have indices in 4 dimensions only, thus the 
algebra is greatly simplified. The D-dimensional indices are introduced by 
the loop integrals, which generate the D-dimensional metric tensor g^^^ so 
that 7^ in D dimensions is given by g^fyp. The basic rules 

g,ug>"' = D , g,,g'^'' = A-D 

g.ug"' = (28) 

are formally similar to the HV ones provided the roles of g^u and ^^^^ are 
exchanged, but, contrary to that case, DRED is known to be algebraically 



inconsistent |2^. Moreover one should also mention that DRED fails to 
reproduce the triangle anomaly, unless further ad-hoc prescriptions are as- 
sumed Finally it is well known that one has to take care of many 



theoretical subtleties regarding the renormalization of operators, living in 
(4 — D) dimensions, when higher order calculations are performed ETf . 



In spite of all these problems, no calculation, with the exception of the 
triangle anomaly, gives so far a wrong result. Concerning our calculation, 
neither inconsistencies of the regularization scheme nor other problems with 
the renormalization procedure in D dimensions show out, once a suitable 
definition of the (4 — D)-dimensional operators is introduced, as explained 
in Sec. |^. 



5 Diagrams and Counter-terms 

In this section the main results of our study are given. Contributions from 
the one and two loop Feynman diagrams in figs. (0)- to the anoma- 
lous dimension matrix are presented in the three considered regularization 
schemes. Tabs. (|l])-(0) contain the results of these diagrams in terms of the 
coefficients of the poles in the number of dimensions, appearing as powers of 
l/e = 2/(4-D). 
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Figure 1: One loop diagrams which generate both the 4-fermion operator 
mixing and the counter-terms to be used in the two loop calculation. 

5.1 One Loop Diagrams 

Let us start with the 0(a°) finite mixing among 4-fermion and magnetic 
operators, already discussed in Sec. ||. The penguin diagrams in fig. (|^) 
could mix 4-fermion operator with magnetic ones. The coupling constant 
is re-absorbed in the definition of the magnetic operators, so this mixing 
appears already at 0{a^). However penguin diagrams in HV and DRED 
have neither pole nor finite parts on the magnetic form factor. Differently 
in NDR, when massive quark fields propagate in the loop, the diagram PI 
induces a finite mixing among Q5, Qq and Q^, Qs- In fact the calculation of 
PI with a 7^ ® 7^_R vertex insertion gives 

2m, (^7^ - 7^^)) (1+ 75). (29) 
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diagram 


M 


7l ® 


7l ® 7mR 






(lA) 


(lA) 


Vi 


2 


1 


1 


Vi 


2 


-4 


-1 




2 


1 


4 




1 


-4/3 






1 


-4/3 


-4/3 



Table 1: Singular parts of the diagrams in figs. (|Tp, with a or 
a 7^ ® 7^i? vertex insertion. V1-V3 results are proportional to the inserted 
4-fermion structures. V4-V5 are proportional to 7^ ® 7^. The multiplicity 
of the diagrams is also reported in the table. Colour factors and ag/'^Ti are 
omitted. 



Selecting the magnetic form factor through the projection 

+ 75)^% ^ ^, + 75)%^ ^ (30) 



the values of the vector Z in eq. (|T2|) can be readily obtained. The accidental 
vanishing of Z in both the HV and DRED schemes makes them to give 
coincident results even at intermediate stages of the calculation. 

Figs. (P and (|]) show all the diagrams required to calculate at the 
leading order the 4-fermion operator mixing matrix 7^ and the mixing of 
the magnetic operators among themselves. The corresponding results can be 
found in tabs. (|I|) and where the pole coefficients are reported. These 
results, which do not depend on the regularization scheme, are established 
since a long time [T^, thus we omit the details on their calculation. Our 
results, as well as the corresponding anomalous dimension sub-matrices in 
eqs. (^31) and (^i]) , coincide with those of refs. |2^ and [Q. We just mention 



that the gauge dependence of the diagrams M1-M5 cancels out when the 
external gluon field renormalization, see eqs. (p7|)-(^), is taken into account. 
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Figure 2: Diagrams responsible for the one loop renormalization of the mag- 
netic operators. 



5.2 One Loop Counter-terms 

Counter-terms to be inserted in the two loop diagrams are obtained from the 
diagrams in figs. Many scheme dependent operators, vanishing in 4 

dimensions, are generated along with the usual counter-terms proportional 
to the operators already present in the 4-dimensional basis. It is well known 



1^, ^ that "effervescent" operators must be considered in the renormal- 
ization procedure, namely they must be properly inserted as counter-terms 
in the two loop diagrams, in order to have the right final result. 

Here we do not repeat the procedure to define these operators in HV 
and NDR, since it is fully explained in refs. [0, |T|, |2D[. Also the DRED 
"effervescent" operators generated by the 4-fermion diagrams in fig. (|lD are 



known from ref. We confirm all previous results and give some details 
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Diagram 




Ml 




M, 


-4 


M, 


-1 


M, 


6 


Ml' 


5 


M, 


3 
2 


5 


-2 



Table 2: Singular parts of the diagrams in fig. @. All the results are 
proportional to the magnetic operators. The factor as/47r is omitted. The 
diagrams M4 and M5 are calculated both in the Feynman and background 
gauges. 



only on our DRED calculation of the "effervescent" operators due to the 
penguin and gluon-photon diagrams in fig. recently presented also in 
ref. [|1T|. 



Contrary to other schemes, in DRED the 7 algebra is performed in 4 
dimensions so that no complicated tensor structure appears and the "effer- 
vescent" operators can be readily defined by inspection as those terms that 
have (4 — D)-dimensional Lorentz indices saturated on the external fields. 
The "effervescent" part of the penguin diagrams in DRED is then given by 

-^^^"^^7.(1- 75). (31) 



This result coincides with the corresponding one of ref. [11 



Concerning the diagrams Cia and Cu in fig. their "effervescent" 
parts are found to be respectively 

~^{q + 2iY^,r''{i - 75) + - {ir - Yi) iur'{i - 75) + . . . 

^{q + 2iriug''P{l - 75) + - ih' - lug'-'il - 75) + . . . (32) 

The dots indicates further terms that can be omitted since they contribute 
only to the "effervescent" part of the two loop diagrams. By summing the 
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Figure 3: Diagrams which generate the counter-terms discussed in the text. 
In DRED they generate "effervescent" contributions that cannot be omitted. 



two terms in eq. (|32|), we obtain a counter-term with a pole coefficient which 
disagrees from the one presented in ref. 0]. We beheve that there is a 



misprint in ref. |]TT|, since we obtain the same results reported there for the 
insertions of this "effervescent" operator in the two loop diagrams. 

It is worthwhile to note that, as shown in ref. [|I^], the insertions of the 
counter-terms generated by the diagrams in fig. (|^) and by the longitudi- 
nal parts of the penguin diagrams vanish, so that one can retain 4-fermion 
operators only as counter-terms in the Hamiltonian. However in DRED the 
diagrams of fig. (H) give "effervescent" contributions that must be taken into 
account. 
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5.3 Two Loop Diagrams 



The relevant two loop diagrams are shown in fig. (^) and the corresponding 
results in the three considered regularization schemes can be found in tabs. 
(|3D-(p!0D, which contain the pole coefficients of the diagrams projected on the 
magnetic form factor. Two tables for each different kind of Dirac structure 
inserted in the upper vertex are presented. The first one contains the pole 
coefficients of the bare diagram (D), the insertion of 4-dimensional counter- 
terms (C) and the insertion of "effervescent" counter-terms (E). The second 
one collects the final results of the renormalized diagrams (D), obtained as 



D = D-C 



2 ' 



(33) 



see eq. (^) below. Results from the insertion of a 7^ ® upper vertex in 
the P-type diagrams of fig. are reported in tabs. Those coming 

from 7^ ® 7^L insertion in F-type diagrams can be found in tabs. 
7l ® 1fj.R 4-fermion vertex, inserted in P- and F-type diagrams, originate 
the results in tabs. and (^)-(|T0|) respectively. Note that the left- 

right operator insertion in the P-type diagrams, tab. (|^), do not vanish just 
because of the massive loop propagators, thus they contribute only when 
q = b is taken in Q5, Qq. 

The calculation of the two loop diagrams result in many tensor struc- 
tures, containing mt, mass and/or external momenta. Further projections 



are required, besides those in eq. 

"^6(1 + 75)7^1*- 



mfe(l + 75)^g^^0 



m,(l + 75)^g'^^0 



mb(l + 75)^^ 



1 
4 



(l + 75)^g'^^0 
(l + 75)|(g^^0 



mfe(l + 75)^7'' - 
^6(1 + 75)^^^^ 
m,(l + 75)^P^ 



1 

4 



4 
1 



(34) 
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Moreover in DRED also the antisymmetric tensor e^jypo- appears, either in 
pure 4-dimensional expressions or in tensors involving (jf^,^. In these cases the 
projections are 

+ 75) W^"""?. ^ 1 
m6(l + 75)7a7p^^'^^^'^P. ^ -7(4-^) 



1 
4 



(l + 75)7a^X'^^%P. ^ -^(4-1^). (35) 

Let us consider now the relations among different schemes. Eqs. (pi]) 
enforce a set of relations among one and two loop diagrams, which can be 
easily obtained by considering how each diagram contributes to the anoma- 
lous dimension matrix (i.e. its Dirac and colour structure). These relations 
are 

A(P2 + P3) = 0, A(F2 + F3) = 
A(P2 + P4) = 0, A(F2 + F4) = 

AP-r = -iPiAP^, AFj = -ipiAPr 

A{Pp + Pp + P^" + Pp + Pp + Pg™") = -TAPr" 

A (Pf' + Pr + Pq™" + Pg™*- ) = -2 AP^ , (36) 

where A indicates the difference between two different regularization schemes. 
Quantities denoted by rrih refer to the results of the left-right P-type diagrams 
with a mass insertion into the loop propagators. They can be found in tab. 
(H), with the exception of P^'', which is equal to 2 in NDR and vanishes 
in the other two schemes. As already noted, HV and DRED calculations 
coincide diagram by diagram, so that eqs. ( p6D are not really interesting. 
However, when comparing HV or DRED with NDR, this check is effective 
and one can verify that all our diagrams satisfy eqs. (B^). 
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Among the other checks passed by our results, one can readily verify 



the cancellation of all the double poles, as an indication that this is a 
leading order calculation, see eqs. (^Tj); 

the usual relation between the double poles of the bare diagram and 
the insertion of the counter-terms, the second being two times larger 
then the first; 

the gauge independence, checked by using the Feynman and the back- 
ground gauges ||26[ . 
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Figure 4: The two loop diagrams relevant for the calculation of the vectors 
A'.s in eq. (|). 
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Table 3: Singular parts of the P diagrams in fig. (^) with a 7^ ® vertex 
insertion. The common double poles and the single poles, calculated in HV, 
NDR and DRED, are presented for the bare diagrams (D), the 4-dimensional 
(C) and the "effervescent" (E) counter-terms. All the results are proportional 
to the magnetic operators. Diagram P4 is calculated both in the Feynman 
and background gauges. 
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Table 4: Singular parts of the renormalized P diagrams in fig. (^) with a 
7i ® 7^L vertex insertion. These results are obtained from tab. (|[) using eq. 

m. 
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Table 5: The same as tab. (|]) for the F diagrams with a 7^ ® vertex 
insertion. 
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Table 6: The same as tab. 
insertion. 



D 



1 

e2 


I HV 
e 


I NDR 
e 


I DRED 
e 


1 


71 


59 


71 


9 


54 


54 


54 


1 


91 


103 


91 


9 


54 


54 


54 


1 


5 


7 


5 


9 


54 


54 


54 


1 


5 


7 


5 


9 


54 


54 


54 




2 


4 


2 




9 


9 


9 



for the F diagrams with a 7^ ® 7^l vertex 
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Table 7: The same as tab. @) for the P diagrams with a 7^ ® 7^j? vertex 
insertion. 
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Table 8: The same as tab. @ for the P diagrams with a 7^ ® 7^/j vertex 
insertion. 
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Table 9: The same as tab. (|^) for the F diagrams with a 7^ 7^/2 vertex 
insertion. 
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Table 10: The same as tab. 
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for the F diagrams with a 7^ ® 7^r vertex 
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6 Anomalous Dimension Matrices 



The anomalous dimension matrix appearing in eq. (^) is defined as 



(37) 



where the diagonal matrices n/ and fig respectively count the number of 
fermion and gluon external fields of the operator they are applied to, {Si)ij = 
SaS^j + SisSsj, and {S2)ij = SisSsj- The anomalous dimensions due to the 
external fields and to the explicit couplings and masses are known to be 



7/ 



I9 



2N 

-2 ( -iV - -Uf 



Air 



7 



^^9 
19 



-6- 



Ar2 



2N 

\3 3 ^ 
N--nf] + 



(38) 



The two values of jg refer to the Feynman and background gauge calculations. 

The operator anomalous dimension 70 is defined in terms of the matrix 
of the renormalization constants as shown in eq. ([T5|). In turn this matrix 
has an expansion, in terms of the renormalized coupling constant as and 
the regularization parameter more involved than in other cases. In fact 
now the renormalization constants include an explicit dependence on the 
subtraction scale /i, starting already at order 0(a°), due to a mismatch in 
the dimension between 4-fermion and magnetic operators. Thus the usual 
expression of the anomalous dimension becomes 



70 



2Z- 



-eas + I3{as)) 



_d_ 

da. 



(39) 



while the multiple expansion of the matrix Z is given by 

1 \ a 



-2e 



^1,2 1 I ^1,2 1 
^1 - + ^2 ~J 

e e'^ 



+ ... 



(40) 
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The coefficients are labeled as Z^'^, where a is the order in as, b is the number 
of loops involved in the calculation and c is the order in the j expansion. 

Using eqs. and (^), we obtain 
Z'/ = 

7o = -2Zl'' -a{^zI'' -Izl'^Z^,'' -Iz'/Z]^ (41) 

In order to have a finite 70 when D tends to 4 dimensions, the first of 
eqs. (|41|) must be satisfied. Then the second equation gives the anomalous 
dimension matrix in terms of the single pole coefficients and finite parts of 
the one and two loop diagrams, see also eq. (|33|). This equation is obtained 
by using the relations 

z''')__ = (zy)__ = o 



'0 IPE ~ V^O /PS 



\ ^ J EP \ ^ J EP \ ^ J EP 



Z'i^^zY = Z'''Z^'' = 0, I, J = 0,1 (42) 



which hold for our MS renormalization constants. P indices get values in the 
set of "physical", i.e. 4- dimensional, operators, while E refers to the "effer- 
vescent" ones. Moreover only matrix elements between "physical" operators 
are retained, the "effervescent" contribution being included in the last two 
terms of 70 in eqs. (^). In fact the summed indices in the products of the Z 
matrices run over the full D-dimensional basis, including the "effervescent" 
operators. 



The expression of 70 in eq. (|41| ) is similar to the next-to-leading order 
formula, see e.g. refs. [^, even if this is a leading order calculation. 
In particular the third term in eq. (PT| ) accounts for the well-known "effer- 
vescent" contribution coming from one loop 4-fermion diagrams. The fourth 
term includes the peculiar operators, present only in the DRED scheme, 
which are generated in 4 — D dimensions by the diagrams in fig. (^. Both 
these terms are actually purely "effervescent" , since they involve finite parts, 
which are not included in the MS renormalization constants as far as matrix 
elements between 4-dimensional operators are concerned. On the contrary. 
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matrix elements connecting "physical" and "effervescent" operators must be 
retained in the matrix of the renormalization constants, even in the MS case. 

In the last part of this section we summarize the results for the anomalous 
dimension matrices in the three considered schemes. The new calculation in 
the DRED scheme gives a matrix which is identical to the HV one. Hence 
the results reported here can also be found in ref. [0]. 

Splitting the anomalous dimension matrix as in eq. we give the 

results for 7^, /^y, /?8 and 777, 737, 788- 

The regularization scheme independent matrix 7^ is given by 



7r 



' N 



6 






6 

__6_ 

N 












2 

"SiV 

22 

~ w 



2n 







2 
3 

22 
3 



6 I 2nf 
TV 3 



2n 
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"37V 
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"sat 

2nt 



3N 

_6 

N 

2nf 





2 
3 

4 
3 

2nf 



-6 



_i2 



3N 



3N 



> N^-1 
2N 



2n 



(43) 



where is the number of colours and = + is the number of active 
flavors. 

The mixing of the magnetic operators O7 and Og among themselves is 
also scheme independent and is given by 



777 
787 



1 



2N 
- 1 

) 

' 2N 



AN 



N' 



(44) 



For the vectors (3, which depend on the regularization scheme, we obtain 
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where nf = nd+ q^^m- The NDR resuh is given by 

/ 



-16 I 12Qu N'^-l 

9 2N ^ Qa 2N 



184 
9 2N 



l&Uf N^-l , 12nf{N^-l) 
'~9 2N ' 2N 



40 



2N 



16n/ Af2-i _ 12nf{N'^-l) _|_ mN{N'^-l) 



2N 



2N 



2N 



(46) 



22N _ m_ 

9 9Af 
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9 m ^ 
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22Nn 



9N 



4iV-f -6n, 



9 9N 



7 Status of the Calculation of the QCD Cor- 
rection to the 6^57 Decay 

Let us start with the original calculations, from ref. to ref. 0. These 
works used the "reduced" basis. As we have shown, this approximation leads 
to scheme dependent results. Apart from this, there is no computational 
error in both NDR and DRED calculations. However refs. |^ did not 
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include the contribution coming from the 0{a^) mixing among Q5, Qq and 
the magnetic operators, which is present in NDR. On the other hand this 
mixing vanishes in DRED. However the authors of ref. [0 overlooked the 
contributions of the "effervescent" counter-terms. These two calculations 
obtained different final results. This difference led the authors of ref. to 
the incorrect conclusion that the DRED scheme fails in this case. Now we 
know that those results can differ (and indeed they do) without implying 
any failure of the regularization scheme. A second, more specific, argument 
in ref. was based on the explicit calculation of the sum of the diagrams 
P2+P3infig. (^). They computed this sum in DRED and in a 4-dimensional 
regularization scheme and obtained again different results. Incidentally eqs. 
( PBD show that the sum P2 + P3 is indeed scheme independent. We have 
checked that the two results of ref. actually coincide, once one includes 
the contribution of the DRED "effervescent" counter-terms. 

Coming to more recent calculations, the full basis and the already men- 
tioned 0{a^) mixing have been taken into account, starting from ref. 0], 
so that now problems with the scheme independence of the Effective Hamil- 
tonian are no more present, as shown in ref. 0. However the three latest 
works on the subject §|, give three different results for some anomalous 
dimension matrix elements. In particular, using our normalization, 

. ref. § gives 757 = -32, = 758 = 10, 768 = 

. ref. i gives 757 = ^, 767 = ^, 758 = -'-f, 768 = 

. ref. i gives 757 = f^, 767 = ^, 758 = -f , 768 = -^ff ; 

when = 3 and Uf = 5. The origins of these differences have been clearly 



summarized in ref. [|lT|. We shortly repeat them here, taking our calculation 
as a reference 

• the calculation of ref. differs from our one because of the values of 
the diagrams P2 and P3 in tab. @; 

• the calculation of ref. [^] differs from our one because the results pre- 
sented there do not include the "effervescent" counter-terms, as also 



stated in ref. 1T2 



31 



Of course refs. |^ and differ one from the other for both the two reasons 
hsted above. 

In our opinion, the results of ref. actually agree with our ones. In 
fact we cannot see any reason not to include the "effervescent" counter- 
terms, which are known to be needed from a long time^, see also eq. (|4l|) . 
Concerning ref. 0, we can just say that our diagrams, including those ones 
responsible for the difference, verify all the checks we have done, including 
those enforced by the scheme independence of the Effective Hamiltonian, eqs. 
(|36D. Up to now we have not been able to find any error in the calculations of 
P2 and P3 with a mass insertion into the loop, which are actually quite 
easy to evaluate. Thus we are confident that our results are correct. 

Anyway these differences still present in the literature are known to have 
a very little impact on the phenomenology of the radiative B decays. We 
plan to present our phenomenological analysis in a forthcoming paper . 
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